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ABSTRACT 



We report on the discovery of > 100 MeV 7 rays from the binary system PSR 
B1259— 63/LS 2883 using the Large Area Telescope (LAT) on board Fermi. The 
system comprises a radio pulsar in orbit around a Be star. We report on LAT 
observations from near apastron to ~ 60 days after the time of periastron, t p , 
on 2010 December 15. No 7-ray emission was detected from this source when it 
was far from periastron. Faint 7-ray emission appeared as the pulsar approached 
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periastron. At ~ t p + 30d, the > 100 MeV 7-ray flux increased over a period of a 
few days to a peak flux 20-30 times that seen during the pre-periastron period, 
but with a softer spectrum. For the following month, it was seen to be variable 
on daily time scales, but remained at ~ 1 — 4 x 10~ 6 cm" 2 s" 1 before starting 
to fade at ~ t p + 57d. The total 7-ray luminosity observed during this period is 
comparable to the spin-down power of the pulsar. Simultaneous radio and X-ray 
observations of the source showed no corresponding dramatic changes in radio 
and X-ray flux between the pre-periastron and post-periastron flares. We discuss 
possible explanations for the observed 7-ray-only flaring of the source. 

Subject headings: pulsars: individual (PSR B1259-63) — (stars:) binaries (PSR 
B1259-63/LS 2883) — gamma rays: stars — X-rays: binaries 



1. Introduction 



The pulsar system PSR B1259— 63 was discovered at Parkes in 1989 and comprises a 
47.76 ms radio pulsar orbiting a massive star (LS 2883) in a highly elliptica l (e ~ 0.87) orbit 
with a period of ~ 3.4 years (I Johnston et al.lll992l ; iNegueruela et al.ll201ll ). Recent optical 
spectroscopy (INegueruela et al.ll201ll ) yields an updated distance estimate to this source of 
2.3 ± 0.4 kpc, in reasonable agreeme nt with the dispersion measure (DM) derived distance 
of 2.7 kpc using the NE2001 model JCordes & Lazioll2002h . so we adopt D = 2.3 kpc. The 
companion shows evidence for an equatorial disk in its optical spectrum, and has generally 
been classified as a Be star ( IJohnston et al.lll994[ ). The pulsar comes within ~ 0.67 AU of its 



comp anion star at periastron, which is roughly the size of the equatorial disk ( 1 Johnston et al. 



19921 ). The orbital plane of the pulsar is believed to be highly inclined with respect to this 
disk and s o the pulsar crosses the disk plane twice each orbit, just before and just after 
periastron (IMelatos et al.lll995l ). Shock interaction between the relativistic pulsar wind and 
the wind and photon field of the Be star is beli eved to give rise to the variable unpulsed 



X-ray emission observed throughout the orbit (jCominsky et al.l Il994l ; IChernyakova et al. 
20091) and the unpulsed r adio and TeV 7 rays observed within a few months of periastron 
( IChernyakova et al.l 120061 ). 



At energies around 1 GeV, the Energetic Gamma-Ray Experiment Telescope (EGRET) 
provided only an upper limit fo r the 1994 periastr on passage (F 7 < 9.4 x 10~ 8 cm -2 s _1 for 
E > 300 MeV, 95% confidence, (ITavani et al.lll996l )). In TeV 7-rays the system was detected 
during the 2004 and 2007 periastron passages and flux variations on daily timescales were 
seen for energies > 0.38 TeV in 2004 (kharonian et alil2005l . l2009h . 
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For the 2010/2011 passage the time of periastron t p was on 2010 December 15. By 
comparison to previous passages, the unpulsed radio and X-ray emission was expected to 
start rising in mid 2010 November peaking around t p — lOd in the pre-periastron phase and 
reaching another peak around t p + 15d in the post-periastron phase. By 2011 April these 
emissions are expected to go back to their levels when the pulsar is far from periastron. 



Abdo et al.l (l2010d ) reported the first discovery of GeV 7-ray emission from this system 



Kone et al. 


(2011 


) and in 


Abdo et al. 


(2011) 



Tarn et al.l (1201 ll ) reported with further details the GeV 7-ray activity from this system. 



We have assembled a multiwavelength campaign to monitor the system in radio, optical, 
X-rays, GeV, and TeV 7-rays during the 2010/2011 periastron passage. Here we describe 
the Fermi-LAT detection of PSR B1259— 63 in the E > 100 MeV range. We also present a 
preliminary analysis of a portion of the radio and X-ray data to determine if there was any 
anomalous multiwavelength behavior compared to previous periastron passages. We have 
analyzed LAT data over the entire time period from the beginning of the Fermi mission 
(2008 August 4; at which time the pulsar was nearing apastron) through periastron up until 
2011 April 22 which is after the passage of the pulsar through the dense equatorial wind of 
the massive star. Full analyses and interpretation of the multiwavelength data are deferred 
to subsequent papers. 



2. Observations and Data Analysis 

Analysis of the Fermi LAT data was performed using the Fermi Science Tools 09-21- 
00 release. The high-quality "diffuse" event class was used together with the P6_v3_diffuse 
instrument response functions. To reject atmospheric 7-rays from the Earth's limb, we se- 
lected events with zenith angle < 100°. We performed standard binned maximum likelihood 
analysis using events in the range 0.1-100 GeV extracted from a 20° x 20° region centered on 
the location of PSR B1259— 63. The model includes diffuse emission components as well as 
7-ray sources within 20° of the source (based on an internal catalog created from 18 months 
of LAT survey data). The Galactic diffuse emission was modeled using the gll_iem_v02 
model and the isotropic component using isotropic_iem_vO2 0. To better constrain the dif- 
fuse model components and the nearby sources, we first generated a model using two years 
of data between 2008 August 4 and 2010 August 4, a period during which the pulsar was 
far away from periastron. We fixed spectral parameters of all the sources between 5° and 



1 



http: / /fermi. gsfc.nasa.gov/ssc/data/ analysis / 
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15° from the source, and left free the normalization factor of all the sources within 5° that 



were flagged as variable source in the 1FGL catalog (lAbdo et al.l l2010af ). Normalizations 
for the diffuse components were left free as well. For this time period, the source was not 
detected with the LAT and we place a 95% upper limit on the photon flux above 100 MeV 
Fioo < 9 x 10" 9 cm" 2 s" 1 assuming a power law spectrum with a photon index T = 2.1. 

The results of this fit were used to constrain the background source model for analyses 
on shorter timescales starting in November 2010. In the source model, the normalization 
of the isotropic component was fixed to the 2-year value, while the normalization for the 
Galactic diffuse component and three variable sources were left free. 

We searched for 7-ray emission from this source on daily and weekly time scales during 
the first disk passage (mid November to mid December 2010). No detection at the level of 
5o" was observed from the source on these time scales. Integrating from t p — 28d (the typical 
start of enhanced X-ray and unpulsed radio flux) to periastron yielded a clear detection 
of excess 7-ray flux from the sourc e with a test statist ic (TS) of ~ 24 which corresponds 



to a detection significance of ~5<r ( iMattox et al.lll996l ). To estimate the duration of this 



enhanced emission and to get the best fit for the spectrum we looked at the cumulative TS 
as a function of time for integrations starting at t p — 28d (Figured]). Inspection of this plot 
reveals that the TS drops monotonically for integrations ending after t p + 18d, so we use this 
as the end of the integration for this initial period of detected 7-ray flux. For the rest of the 
paper we will refer to this period as the "brightening" . 

During the brightening, the detected 7-ray signal was in the energy range 0.1-1 GeV and 
no significant emission was detected above 1 GeV. The spectrum in the energy range 0.1-1 
GeV is best described by a simple power law with photon index T = 2.4 ± 0.2 sta t ± 0.5 sys 
with average photon and energy fluxes of F = (2.5 ± 0.8 stat ± 0.8 sys ) x 10~ 7 cm -2 s -1 and 
F = (0.9 ± 0.3 stat ± 0.4 sys ) x 10" 10 erg cm" 2 s" 1 respectively. Because of the low signal-to- 
noise ratio during this period spectral fits to an exponentially cutoff power law were not 
constraining. This period is shown as the shaded region on the top panel of Figure [21 which 
shows the 7-ray flux in weekly time bins in the period t p — 131d to t p + 128d. 

At about t p + 30d the source brightened rapidly, reaching a flux that was ~ 10 times 
higher than the integrated flux measured during the first disk passage (Figure [2]). This 
flare lasted for about 7 weeks. The spectrum during this period (t p + 30d to t p + 79d) is 
best described by a power law with exponential cutoff. The best-fit result is obtained for 
a photon index T = (1.4 ±0.6 sta t ± 0.2 sys ) and a cutoff energy E c = (0.3 ±0.1 stat ± 0.1 sys ) 
GeV. The average photon and energy fluxes above 100 MeV during this period are F 10 o = 
(1.3±0.1 stat ±0.3 sys )x 10" 6 cm" 2 s" 1 and G 100 = (4.4±0.3 stat ±0.7 sys ) x 10" 10 erg cm" 2 s" 1 
respectively. 
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In addition to the significant difference in flux and spectral shape between the brightening 
and the flare periods, the weekly time bins show clear evidence for a change in the photon 
index for a power-law fit on these time scales (bottom panel of Figure H]). The phonon index 
softens from 2-2.5 during the brightening to a value of 3.5 around the peak of the flare. After 
that, the index hardens over the rest of the flare period to its values during the brightening. 



The top panel of Figure E] shows 7-ray flux as a function of time in 1-day time bins 
during the flare. The daily source flux du ring this per i od is y ariabile at the 99.9% confidence 
level according to the method outlined in lAbdo et al.l ( I2010af ). During this strong- variability 
flaring period the source flux varied by a factor of 2-3 on daily time scales. 

The previous upper limits from EGRET are fully c onsistent with the flux observed by 
Fermi during the same orbital phase (ITavani et al.lll996l ). At the time of the bright emission 
seen by Fermi, which was well above the sensitivity level of EGRET, EGRET was pointed 
elsewhere. 



3. Timing Analysis 

To search for evidence of 7-ray pulsations from PSR B1259— 63, we have constructed 
a radio ephemeris using observations from the 64-m Parkes telescope. Using Tempo2 



( jHobbs et al.l 120061 ) . we fitted a timing model to 45 TOAs covering the range 2007 Oct 
6 through 2011 January 16. Using this ephemeris, we folded LAT photons in the interval 
2008 August 4 through 2010 August 4, a total of 24 months of observation, all well away from 
the flaring region near periastron where non-pulsed 7-ray emission is observed. We found 
no statistically significant indication of a pulsed 7-ray signal from this source. The pulsed 
flux upper limit depends on the unknown pulse shape and assumed spectrum. We therefore 
use the continuum upper limit on the energy flux above 100 MeV of G100 < 1.0 x 10 -11 erg 
cm -2 s _1 in our comparisons with the rest of the 7-ray pulsar population. 

Comparing this pulsar to the rest of the LAT-detected pulsars, we find that most de- 
tectability metrics predict that this should be a 7-ray pulsar. Although the characteristic 
age of 333 kyr is fairly large, the spin period is short for a middle-aged pulsar and thus at a 
distance of 2.3 kpc the E of 8.2 xlO 35 erg s _1 and magnetic field at the light cy linder of 2.9 



xlO 4 G are well within the range where 7-ray pulsations are typically detected (lAbdo et al. 



2010bl ). If we assume that the beaming factor fa is 1, then the 7-ray efficiency of the pulsar 
is less than 0.7%, which is considerably lower than the 10% efficiencies that are typical in 
this E range. Determining whether this pulsar is intrinsically under-luminous in the 7-ray 
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band or if the low 7-ray luminosity is simply a geometric effect will require detailed modeling 
that includes geometrical information from radio polarization measurements. 

We have also searched for 7-ray pulsations during the brightening and flare periods where 
continuum 7-ray emission was detected. No pulsations were detected in these intervals. This 
is consistent with this 7-ray emission originating from the intrabinary shock, which is well 
outside the light cylinder of the pulsar and thus is not expected to be modulated at the spin 
period. 



4. Radio and X-ray monitoring 



Pulsed emission was monitored at Parkes to look for changes in the DM and rotation 
measure (RM) and determine the duration of eclipse of the pulsed signal. Pulses disappeared 
on t p — 16d and reappeared on t p + 15. In the ~ 2 weeks leading up to the disappearance of 
the pulse, significant changes in the DM were observed. 

The PSR B1259-63 system was monitored at frequencies between 1.1 and 10 GHz using 
the ATCA array. Twelve observations spanning t p — 31d to t p + 55d were collected. Unpulsed 
transient radio emission was detected throu ghout the periastron passage with a behavior 
similar to that seen in previous observations (I Johnston et al.ll2005l ). as shown in Figure HI 



X-ray observations of PSR B1259-63 during this passage demonstrated the repeatability 
of the 1-10 keV light curve as shown in FigureHl As with the periastron passage of 2004, Swift 
observed a rapid X-ray brightening starting at ~ t p — 25d. These observations confirmed the 
spectral hardening precedi ng the pre-periastron flux rise. Similarly to previous periastron 
passages (see for example IChernyakova et al.l ( 120091 ) and references therein), observations 
with Swift , Suzaku and XMM-Newton showed a rise of X-ray flux after periastron. 



5. Discussion 



Emission from the PSR B1259-63 system is produced in the interaction of the pulsar 
wind w ith the stellar wind of the companion star. Observa t ions in radio, X-ray and TeV 7-ray 



bands (Johnston et al 



1992 



19961 ; IChernyakova et al. 



2006 



Aharonian et al.ll2005l . l2009t iKawachi et al.l 12004c iTavani et al. 



20091 ) revealed a characteristic variability of this emission 



during the periods of periastron passage. Detection of the 0.1-10 GeV band 7-ray emission 
around periastron was not unexpected. However, Fermi observations reveal puzzling behav- 
ior of the source, which was not predicted in any model of 7-ray emission from this system. 
An unexpected strong flare, visible only in the GeV band was observed some 30 days after 
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the periastron passage and after the neutron star passage of the dense equatorial wind of 
the massive star. 

During this flare the source was characterized by an extremely high efficiency of con- 
version of pulsar spin-down power into 7-rays. The highest day-average flux was F 10 o ~ 
3.5 x 10~ 6 cm -2 s _1 with a spectral index of T ~ 3.0. This corresponds to an isotropic 7-ray 
luminosity of ~ 8 x 10 35 (-D/2.3 kpc) 2 erg s , nearly equalling the estimated total pulsar 
spin down luminosity Lsd — 8.3 x 10 35 erg s _1 ( 1 Johnston et al.lll992l ). This is illustrated in 
Figure where the horizontal red line shows the flux which would be produced when 100% 
of the spin-down power is converted into radiation emitted within one decade of energy, not 
taking into account possible beaming effects. 

Broadband spectra of emission around periastron are shown in Figure Strong in- 
creases in GeV flux and changes in 7-ray spectrum during the flare were not accompanied 
by noticeable spectral variations in the X-ray band. 

Several possible mechanisms of production of 0.1-10 GeV 7-ray emission from the sys- 
tem were previously discussed: synchrotron, inverse Compton (IC), Bremsstrahlung, or 



pion decay emission (|Tavani fc Aronsl 119971 ; iKawachi et al. 



20041 : IChernvakova et all 12006 



100 TeV produce synchrotron emis- 
2 



Khangulyan et al.l 120071 ). Electrons with energies E e 

sion in the energy range E 1 ~ 10 9 [B/l G] [E e /10 u eV] z eV. Alternatively, electrons with 
energies E e ~ 1 — 10 GeV could produce 7 quanta with energies E 1 ~ 10 8 [E e /1 GeV] 2 eV 
via IC scattering of Be star photons. Bremsstrahlung emission in the GeV band could be 
produced by the GeV electrons. Finally, the dense equatorial stellar wind could provide a 
sufficiently dense target for proton-proton interactions followed by decays of neutral pions 
into photons. 

Figure E] shows example model fits to the persistent emission data. The model shown in 
the upper panel ass umes that high-energy part i cles es cape with the speed of th e stella r wind, 
as in the model of IChernvakova Sz Illarionovl (1l999h and IChernvakova et al.l ( 120061 ). Slow 
escape of the high-energy particles leaves enough time for the efficient cooling of electrons 
via IC and/or Bremsstrahlung / Coulomb loss mechanisms. In the lower p anel, high-energy 
partic les are assumed to escape with the speed 10 10 cm s _1 , as in the model of iTavani fc Arons 
( 119971 ). In this case only synchrotron cooling is efficient. The code used for the calculations 
is described in IZdziarski et al.l (120101 ) . 



In general, the flare could be explained either by anisotropy of the 7-ray emission, or 
by an abrupt change of physical conditions in the emission region or by the appearance of a 
new emission component. 



Several possible sources of anisotropy are present in the system: relativistic beaming 
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of the 7-ray emitting outflow from the system (IBogovalov et al.l 120081 ; iDubus et al.ll2010l ). 
anisotropy of pulsar wind, or anisotropy of radiation field of the massive star. The model 
for the flare spectrum (cyan data points) shown in the lower panel of Fig. assumes a 
particular type of anisotropy which could appear at the high-energy end of the synchrotron 
spectrum, at the energies at which the synchrotron cooling distance is comparable to the 
gyroradius. In such a situation the electron distribution could not be isotropized within 
the synchrotron cooling time scale. The assumpiton that highest energy electrons with 
anisotropic initial velocity distribution cool before being isotropized, results in the increase 
of apparent luminosity by a factor 47r/f2 ~ 10 where Q is the solid angle into which most of 
the highest energy synchrotron power is emitted. 

Another possibility to explain the flare is considered in the model shown in the upper 
panel of Fig. It iassumes a local increase of the density of stellar wind by a factor of ~ 10 
which results in the increase of the Bremsstrahlung component of emission spectrum. 

Clarification of the physical mechanism of the puzzling flare discovered by Fermi requires 
more complete view of the properties of the flare, including information on system behaviour 
in optical and TeV bands and on orbit-to-orbit variations of the GeV flaring pattern. 
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Fig. 1. — Cumulative TS from integrations starting on t p — 28d (see text). The shaded area 
marks the 10 days during which the LAT was in a modified rocking mode giving ~ 30% 
higher exposure on the source. The arrow marks the time we adopt as the end of the first 
emission period. 
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Fig. 2. — Gamma-ray flux and photon index of PSR B1259-63 in weekly time bins between 
t — 131d to t+ 128d. Upper panel: > 100 MeV flux, 2er upper limits are drawn for points with 
TS < 5. Bottom panel: variations of spectral index of a power law spectrum. The shaded 
area shows the brightening period. Dashed line marks the time of periastron. Dashed- 
dotted lines marks the orbital phase during which EGRET observed this source in 1994 
( iTavani et al.lll996l ). 
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Fig. 3. — Gamma-ray flux and photon index of PSR B1259-63 in daily time bins during the 
flare. Upper panel: > 100 MeV flux, bottom panel: spectral index of a power law spectrum. 
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Fig. 4. — Light curves of PSR B1259-63 around periastron. Panel a: HESS 2004 and 
2007 periastron passages (lAharonian et al.ll2005l ). Panel b: Fermi -12010 periastron pas- 
sage. Panel c: X-ray flux es from three periastron passages in units of 10 -11 erg cm -2 s" 1 
( jChernyakova et all 120091 ) . Panel d: R adio (2.4 GHz) flux de nsities measured at ATCA for 
the 2010 and 1997 periastron passages (I Johnston et al.lll999l ). 



- 18 - 




-6-5-4-3-2-1 1 2 3 4 5 6 7 8 9 10 11 12 13 

log E [eV] 



Fig. 5. — Spectral energy distribution of PSR B1259-63 around periastron. Blue and cyan 
points represent the measurements of the spectra in the pre- and post-periastron periods (la- 
belled 1 and 2 respectively) by the LAT, Swift-XBT in X-rays and ATCA in radio. Thin solid, 
dotted and dashed lines show synchrotron, inverse Compton and Bremss trahlung compo 



nents correspondingly. Green points show HESS measurements from 2004 ( lAharonian et al. 



20051 ). Light grey curves show the models of pre-periastron emission, dark grey curves show 



the models of the flare. The horizontal red mark shows the flux which would be produced if 
100% of the pulsar spin-down power were converted into electromagnetic radiation. In the 
upper panel the high-energy particles are assumed to escape from the system with the speed 
of the stellar wind, while in the lower panel, high-energy part icles are assumed to escape 
with the speed c/3, as in the model of iTavani fc Aronsl ( 119971 1. see text for details. Idata 
points will be made available through https://www-glast.stanford.edu/cgi-bin/pubpub 



